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ABSTRACT 

Context. The origin of the so called 'detached shells' around AGB stars is not fully understood, but two common hypotheses state that 
these shells form either through the interaction of distinct wind phases or an eruptive mass loss associated with a He-shell flash. We 
present a model of the formation of detached shells around thermal pulse asymptotic giant branch (TP-AGB) stars, based on detailed 
modelling of mass loss and stellar evolution, leading to a combination of eruptive mass loss and wind interaction. 
Aims. The purpose of this paper is first of all to connect stellar evolution with wind and mass loss evolution and demonstrate its 
_C consistency with observations, but also to show how thin detached shells around TP-AGB stars can be formed. Previous attempts to 

Ck link mass loss evolution with the formation of detached shells were based on approximate prescriptions for the mass loss and have 

not included detailed modelling of the wind formation as we do here. 

Methods. Using stellar parameters sampled from an evolutionary track for a 2M Q star, we have computed the time evolution of 
the atmospheric layers and wind acceleration region during a typical thermal pulse with detailed radiation hydrodynamical models 
including dust formation. Based on these results, we simulate the subsequent circumstellar envelope (CSE) evolution using a spherical 
hydrodynamic model. 

Results. We find that existing simple mass loss prescriptions all suggest different mass loss evolutions and that they differ from our 
detailed wind modelling. The most important factor for the formation of a detached shell is the wind velocity evolution which has a 
strong impact on the wind interaction and the resulting pile-up of matter. Our CSE model shows that a thin shell structure may be 
' formed as a consequence of a rather short phase of intense mass loss in combination with a significant variation in the wind velocity, 

as obtained by our wind models. This situation can only be obtained for a limited range of amplitudes for the piston boundary used in 
, the dynamic atmosphere models. 

■ Conclusions. The combined mass loss eruption and wind interaction scenario for the formation of detached shells around AGB stars 

(suggested by previous work) is confirmed by the present modelling. Changes in mass loss rate and wind velocity due to a He-shell 
flash are adequate for creating distinct wind phases and a 'snow plow effect' that is necessary to form a geometrically thin detached 
shell. The derived properties of the shell (i.e. radius, thickness and density) are more or less consistent with existing observational 
constraints. 
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5r . 1 ■ Introduction A more common, older hypothesis for how these structures 

are formed, states that they arise from the interaction between 

The origin of the so called 'detached shells' around TP- a fast and a slow wind (two-wind interaction) from the cen- 

AGB stars in a late evolutionary stage is not fully under- tral star. As the fast wind propagates outwards, it will sweep 

stood. As a phenomenon, the se structures have b een known up material ejected through the slow wind and create a shell- 

for about 20 years (see e.g. IQlofsson et al. 1987I ). The idea like structure. This shell then contains matter ejected by the fast 

that they were connected with He-shell flashes soon emerged, wind, together with swept up matter ejected by the slow wind, 

but t he exact formation m echanism remained a matter of de- In order to work, this hypothesis needs a sufficiently high mass 

bate dOlofsson et al. 19901 ). In fact, it is still not clear that de- loss rate and wind velocity associated with the fast wind, as 

tached shells are only of internal origin. Wareing et al. (2006) well as a non-negligable pre-flash wind (for more details, see 

have recently shown that shell-like structures like these can ISteffen et al. 1998I IQlofsson et al. 2000I ISchoier et al. 20031 and 

emerge from the interaction of an AGB wind with the interstellar references therein). Kwok et al. (1978) presented an analytic 

medium (ISM) in which the star is moving. The shell structure wind interaction model in which a fast wind is running into a 

could then be explained as a bow shock into the ISM. slower wind - each having a constant mass-loss rate and con- 
stant velocity. Hence, the distance from the central star to the in- 

teraction region is expected to increase linearly over time. This 

Send offprint requests to: Lars Mattsson simple picture has been used to link observed mass loss proper- 
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Another early hypothesis (Zuckerman 1993) for the forma- 
tion of thin shell structures suggests that a short event of high 
mass loss (with a non-evolving wind velocity), which we re- 
fer to as a "mass loss eruption", would simply translate into 
a correspondingly narrow, high-density shell moving through 
the circumstellar envelope. This, however, was shown to be 
an insufficient mechanism for the formation of detached shells 
dSteffen et al. 1 9971 [Stiffen et al. 1998b . But the same mass loss 
evolution could still produce the detached shell phenomenon if 
the wind before the flash is slower than the wind during the erup- 
tion, due to the sweep-up of matter that occurs as the matter re- 
leased by the mass loss eruption propagates through the CSE cre- 
ated by the previous mass loss. This mass loss eruption scenario 
is most likely tied to a He-shell flash (thermal pulse) as sug- 
gested by Olofsson et al. (1990) and Vassiliadis & Wood (1993). 
Under the assumption that the slow (pre-fiash) stellar wind rep- 
resents an ordinary state of mass loss for a TP-AGB star there 
is still the need of a period of relatively high mass loss and an 
enhanced outflow velocity, also in an interacting wind scenario. 

Previous numerical models of the connection between a He- 
shell flash and the formation of a detached shell have pro- 
vided some support for the wind-interaction hypothesis, but 
as mentioned above, only as a partial explanation. Steffen & 
Schonberner (2000), from hereon abbrevated SS2000, confirmed 
this picture of a detached shell forming from a mass loss eruption 
and subsequent wind interaction qualitatively, using a spherical 
hydrodynamic model, altough the mass loss evolution is pre- 
scribed, not modelled, i.e. a more or less ad hoc assumption 
about the inner boundary. 

A rather detailed 3D model of the hydrodynamical features 
of a detached shell was presented by Myasnikov et al. (2000), 
where the gasdynamic features of detached shells around a car- 
bon rich TP-AGB star with variable mass loss rate were inves- 
tigated numerically. They find that these shells are unstable to 
turbulence and that Rayleigh-Taylor instabilities develop in the 
radiative case. Furthermore they find that 2D calculations show 
quite different features compared with the 3D case and that the 
structures of the predicted shells are dependent on the initial ba- 
sic wind properties. The latter is a very important finding. In 
order to confirm the two-wind interaction hypothesis on a more 
detailed level, it is thus cruical to connect the evolution of the 
stellar parameters (M, L*, T^) with the mass loss and wind evo- 
lution during the He-flash event. 

To understand the mass loss evolution of TP-AGB stars it 
is important to know how stellar evolution affects the physics 
of the atmosphere and consequently the formation of a stellar 
wind. This can be done by feeding the stellar parameters result- 
ing from a stellar evolution model into a dynamic stellar atmo- 
sphere model. The models for the atmospheres and winds used 
in the present paper are computed using frequency-dependent 
radiative transfer for the gas and dust, including detailed micro- 
physics of the dust grains and their formation (described in 
Hofn er et al. 20031 and lAndersen et al. 20031 in combination 
with time-dependent hydrodynamics. It has been shown that syn- 
thetic opacity sampling spectra for carbon-rich AGB stars based 
on these models reproduce the observed spectral energy distri- 
butions of these stars quite well (Gautschy-Loidl et al. 2004). 
Furthermore, for the first time, these models are also able to si- 
multaneously reproduce the observed time-dependent behaviour 
of fundamental, first and second overtone rotation-vibration 
lines of CO, which are formed in the outflow, wind accelera- 
tion region and atmosphere respectively (Nowotny et al. 2005a, 
2005b). The mass loss properties, however, are more difficult 
to compare with observations. Derivations of 'observed' rates of 



mass loss (e.g. Sch oier et al. 20 05 1 rely on stellar parameters and 
absolute distances that can be uncertain in many cases, which 
means that observations can only provide limited constraints on 
mass loss rates and dust-to-gas ratios. Modelling detached shells 
might therefore be an independent test of the mass loss proper- 
ties of our dynamic atmosphere model. 

In this paper we combine this radiation-hydrodynamics 
(RHD) model of a pulsating TP-AGB atmosphere and wind ac- 
celeration region with a detailed stellar evolution model and a 
simple spherical hydrodynamic model to follow the dynamical 
evolution of the CSE on long time scales and study the forma- 
tion of a geometrically thin detached shell around the star. We 
are thus dealing with three different scales: the internal evolu- 
tion, the atmosphere and the effects on the CSE. As this paper 
shows, our theoretically expected wind properties are such that 
a thin detached shell is formed as a natural consequence of a 
He-shell flash. 



2. Theory and Method 

2. 1 . General Description 

On a descriptive level, the basic idea behind the modelling is 
to connect models of three different regimes (see Table [TJ. The 
starting point is a pre-computed model of the internal evolution 
of a star with a main sequence mass of 2M Q and an initial metal- 
licity of Z = 0.01. This track is not chosen for any particular 
reason, other than that it is consistent with the type of carbon 
stars with known detached shells. The resulting stellar parame- 
ters (M, L*, r e ff, C/O) are sampled with suitable time intervals 
around a thermal pulse. These parameter values are then fed into 
an RHD atmosphere model that provides us with the mass loss 
properties. Knowing these properties we can construct a time 
dependent inner boundary condition for a hydrodynamic model 
of the larger-scale circumstellar evolution, i.e. the formation of 
a detached shell from the interaction of different wind phases. 
The reason why we split up the modelling in this way, is the 
very different time and length scales of the relevant physical 
processes involved. An 'all-in-one' model would demand unrea- 
sonable computing power in order to resolve the smallest lenght 
scales and the shortest time scales. In this section we describe 
the models and how the computations are carried out. 



2.2. Dynamic Atmospheres 

The ID (spherically symmetric) RHD code predicts the mass- 
loss rates of TP-AGB stars by detailed modelling of the atmo- 
sphere and the circumstellar environment around these pulsat- 
ing long-period variable stars. The model includes frequency- 
dependent radiative transfer and dust formation, i.e. we are 
solving the coupled system of frequency-dependent radia- 
tion hydrodynamics and time-dependent dust formation (cf. 
IHofner et al. 20 03 ) employing an implicit numerical method and 
an adaptive grid. 

In our model the stellar atmosphere and circumstellar enve- 
lope are described in terms of conservation laws for the gas, the 
dust and the radiation field, i.e. we solve the following set of 
coupled, nonlinear partial differential equations (PDEs): 

- The three equations describing conservation of mass, mo- 
mentum and energy for the gas. 

- The 0th and 1st moment equation of radiative transfer. 

- The four moment equations of dust formation (cf. Section 

E3. 
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- The Poisson equation (self gravity). 

To this system of nonlinear partial differential equa- 
tions we add a "grid equation" which determines the loca- 
tions of the grid points according to accuracy considerations 
( Dorfi & Drury 19 87) and an equation keeping track of the con- 
densable amount of carbon, leaving us with a total of 12 par- 
tial differential equations (PDEs) to solve. This system of PDEs 
is then solved implicitly using a Newton-Raphson scheme. All 
equations are discretised in a volume-integrated conservation 
form on a staggered mesh. The spatial discretisation of the ad- 
vection term is a monotonic second-order advection scheme 
( Ivan Leer 1977K The same order of numerical precision is used 
for all PDEs. Details of the numerical method are discussed by 
Dorfi & Feuchtinger (1995) and in several previous papers about 
dust-driven wind models (cf. Ho fner et al. 1 995 and references 
therein). 

The dynamical calculations are carried out as follows. All 
wind models are started from hydrostatic dust-free initial mod- 
els where the outer boundary is located close to the photosphere 
(about 1.3 — 1.6/?*). In these initial models, the fundamental 
stellar parameters (M, L+, T e ff, C/O) can be chosen freely as ini- 
tal conditions for the models. In the present case we pick the sets 
of stellar parameters from an evolutionary track as described in 
Section [2~T| When the dust equations are switched on, dust for- 
mation starts and creates an outward motion of the dust and the 
gas. In the first computational phase the expansion of the atmo- 
spheric layers is followed by the grid to about 20-30 R+ (usually 
~ 10 15 cm). At this radius the location of the outer boundary R oul 
is fixed allowing outflow. The outflow model then evolves for 
typically more than 200 years, but the precise time is dependent 
on the pulsation period. The amount of mass between the inner 
boundary R[ n and the outer boundary R out is roug hly 10" 2 M Q , 
which is sufficient to avoid any significant mass depletion within 
the time series considered here. 

The model is self-consistent as far as the physics of the stel- 
lar atmosphere is concerned, but it does not include a physi- 
cal model for the pulsation mechanism. The pulsation is mod- 
elled phenomenologically using a 'piston boundary condition' 
(Bowen 1988) located at R^ n w 0.9R+. To restrict the number of 
free parameters, we employ an empirical period-luminosity rela- 
tion {Feast et a l 19891 1 and keep the period locked to the luminos- 
ity for all models. The piston amplitude, however, is essentially 
a free parameter that we cannot do away with as easily. For this 
reason we must consider the effects of picking different ampli- 
tudes, which is a complication that we will get back to in Section 
14.11 The motion of the piston is accompanied by luminosity vari- 
ations, since the flux through the lower boundary is kept fixed (as 
described in previous papers, e.g., Ho fner & Dorfi 1997b . Note, 
however, that assuming no artificial scaling of the short-term lu- 
minosity variations due to the pulsations, these are more than a 
factor of two weaker than the long-term luminosity variation due 
to the flash. 

2.3. Dust Formation and Mass Loss 

Mass loss from TP-AGB carbon stars is mainly due to momen- 
tum transfer from radiation to atmospheric dust grains. This is 
often referred to as a dust driven wind. As such, a wind model 
for TP-AGB stars is highly dependent on how dust formation is 
treated. Hence, we will give a somewhat more detailed descrip- 
tion of the dust model. 

The models presented here include a time-dependent de- 
scription of dust grain growth and evaporation using the moment 



method by Gail & Sedlmayr (1988) and Gauger et al. (1990). 
The dust component is described in terms of moments Kj of the 
grain size distribution function weighted with a power j of the 
grain radius. The zeroth moment, Ko, is the total number density 
of grains (simply the integral of the size distribution function 
over all grain sizes), while K3 is proportional to the average vol- 
ume of the grains. The equations, which determine the evolution 
of the dust component, are 

^+V-(tfov) = 3, (1) 

dK i j 1 ,73 

-^l + V-{K jV )= J --K hl+ NfZ, 1 < 7 < 3, (2) 

where v is the flow velocity of the dust component, 3 is the net 
grain formation rate per volume, 1 /r is the net growth-rate of the 
dust grains, and is the lower size limit of grains contributing 
to*,. 

In order to calculate how much of the radiative energy and 
momentum is transferred to dust grains, we need to know the 
frequency-dependent opacity of these grains. This can be ex- 
pressed in terms of the extinction efficiency Q ext , which is the 
ratio of the extinction cross section to the geometrical cross sec- 
tion of the grains. Since the grains are usually small compared 
to the wavelengths, Q ext becomes a simple function of the grain 
radius. The wavelength and grain size dependence of the opac- 
ity can thus be separated into two independent factors, which 
greatly simplifies the calculations. The models in this paper are 
calculated using the data of Rouleau & Martin (1991) for Q ext 
(see lAndersen et al. 2003l and lHofner et al. 2 003 for further dis- 
cussion). The intrinsic dust density used in the model is set to 
Pd = 1.85 g cirT 1 , which matches the material in Rouleau & 
Martin (1991). 

We assume spherical dust grains consisting of amorphous 
carbon only, since the star we are modelling is a carbon star (see 
lAndersen et al. 20031 for details about these assumptions). The 
nucleation, growth and evaporation of grains is assumed to pro- 
ceed by reactions involving C, C2, C2H and C2H2. In this model 
of grain growth, a so called sticking coefficient (sometimes re- 
ferred to as the reaction efficiency factor) is used, that enters into 
the net growth rate of the dust grains. This parameter, as, is not 
definitely known unless we know the exact sequence of chemical 
reactions responsible for the dust formation. However, Gail & 
Sedlmayr (1988) argued that the sticking coefficient must be on 
the order of unity, mainly because it is expected that neutral radi- 
cal reactions play a major role in the formation of carbon grains. 
It should be pointed out as well, that with as = 1, our models 
nicely reproduce the expected mass loss properties of TP-AGB 
stars dSchoier et al. 20011 ISchoier et al. 20051 ). See Section |4~T1 
for some further details. 

Dust grains in a stellar atmosphere influence both its en- 
ergy and momentum balance. For simplicity we assume com- 
plete momentum and position coupling of gas and dust, i.e. the 
momentum gained by the dust from the radiation field is di- 
rectly transferred to the gas and there is no drift between dust 
and gas. However, this strong coupling between the dust and 
the gas phase is not obvious. In a previous attempt to relax this 
phase coupling approximation, Sandin & Hofner (2003, 2004) 
found that the effects of decoupling of the phases might be quite 
significant. The most striking feature is that the dust formation 
may increase significantly, but this does not seem to necessarily 
increase the predicted mass loss rates for a given set of stellar 
parameters. However, these models are not directly comparable 
with observations, since one essential component - frequency 
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Table 1. Schematic overview of the model. Radial and time scales given are approximate figures. 



Model 


Input 


Output 


Radial scale 


Time scale 


Comment 


Stellar evolution 


Mzams, Z 


M„L„r eff , c/o 


R <R* 


~ 10 4 yrs 


Pre-computed 


Atmosphere / wind 


M*,L*,r eff , 
C/0,P,Av p 


M, v wind (r) 


R+ < R < 25R-,, 


~ 10 2 yrs 


Full RHD + dust, 
"snap shots" of evolution 


CSE evolution 


M, v wind 


PenvM, V cnv (» 


25/?* <R,< 10 4 fl* 


~ 10 4 yrs 


Gas dynamics, 

dust dynamics not considered 



dependent radiative transfer - is not included. Comparison with 
observations (|Gautschy-Loidl et al . 2004 Nowotny et al. 2005a 



|Nowotny et al. 2005b) > seems to indicate that in general the in 
elusion of frequency-dependent radiative transfer is crucial for 
a realistic description of the pulsating amosphere, and therefore 
for the conditions in the wind formation region. A coupled so- 
lution of the detailed frequency-dependent equations of RHD, 
including dust formation and drift would lead to a prohibitive 
computational effort and is not a realistic option in the present 
context. Thus, we chose to give frequency-dependent radiative 
transfer priority over relaxing the phase coupling approximation, 
even if this approximation cannot hold for very low gas densi- 
ties, i.e. at large distances from the star or for models with low 
mass loss rates. Furthermore, since transfer of internal energy 
between gas and dust is negligible compared to the interaction 
of each component with the radiative field (Gauger et al. 1990), 
we assume radiative equilibrium for the dust. This allows us to 
estimate the grain temperatures from the radiation temperature 
as we know the dust opacities. 

From the models described in this section and in !2.2l we ob- 
tain the density and the wind velocity, both as functions of radius 
and time. The mass loss rate is then given by 



M = 47rRlpo(t)v (t), 



(3) 



where po and Vo are the density and wind velocity at Rq, i.e. 
the radius where the outer boundary is fixed, which is usually 
around 10 15 cm and enough to ensure that the wind has reached 
the terminal velocity. Since practically all momentum transfer 
from radiation to matter is due to the interaction with dust, we 
expect the wind velocity to be correlated with the degree of dust 
condensation. Thus, Eqn. ® leads to the expecation that under 
most circumstances a larger proportion of dust will correspond 
to a higher mass loss rate. 

2.4. Stellar Evolution 

To model how the atmosphere evolves on longer time scales (sig- 
nificantly longer than the pulsation period) due to changes in the 
internal constitution and nuclear processes of the star, we have 
pre-computed a stellar evolutionary track (see Fig. [T}. The stel- 
lar parameters are then sampled at critical or characteristic points 
(tick-marked in the lower part of each panel in Fig. [TJ around a 
thermal pulse in the evolutionary track and used as input to the 
atmosphere model. In total we pick 21 points along the track that 
we know from experience will capture the important features. 
By computing a set of models in this way for different points 
in the time evolution, we can obtain the change of atmospheric 



properties and wind characteristics (i.e. mass loss rate and wind 
velocity) during the flash event. 

The calculations have been done with the one-dimensional 
stellar evolution code EVOL (Herwig 2000, 2004 and refer- 
ences therein). We use the EVOL implementation of exponen- 
tial, time and depth dependent overshooting, here with an ef- 
ficiency / = 0.016 at the bottom of the convective envelope, 
but no overshooting at the bottom of the He-shell flash con- 
vection zone. EVOL features updated micro-physics input, in- 
cluding a nuclear reaction network. The Schwarzschild criterion 
for convection is used and the mixing length parameter is set to 
ffMLT = 1 -7. Mass loss during AGB evolution is included as pre- 
scribed by Blocker 1995 Specifically we use Eqn. (11) from that 
work in which T{M) is determined according to his Eqn. (13) 
and (14). For the mass in Eqn. (13) the intial ZAMS mass is as- 
sumed as an approximation of the actual mass at which the oscil- 
lation period is Pq = lOOd. Mr in Eqn. (1 1) is the Reimers mass 
loss rate (Re imers 19751 1 as given by Eqn. (1) in lBlockeF 1995 
The scaling factor in that equation is set to 77R = 0.1, and gives 
an overall evolution which is consistent with many obserational 
contraints, e.g., the number of thermal pulses expected after 
transformation into a carbon star. 

For this detailed study we have picked the last but 
one (the 17th out of 18) thermal pulse of sequence ET13 



(Herwig, Austin & Lattanzio 2006b. This sequence has the met- 



alicity Z = 0.01 and an initial mass of 2M Q . At the time of the 
selected thermal pulse the total mass is 1.22M and the core 
mass is O.59M . The third dredge-up is well developed after the 
selected thermal pulse, with an efficiency of A = 0.154, after 
A = 0.372 at the previous pulse. Dredge-up efficiency decreases 
at the end of the TP-AGB evolution when the envelope mass de- 
creases. 

To minimise the number of atmosphere models needed to 
obtain the time evolution, we assume a constant stellar mass 
(M* = 1.2M ) during the thermal pulse. This is a reasonable 
simplification in this context even if the stellar mass will de- 
crease notably on longer time scales due to the steady mass loss 
before the first flash and in between subsequent flashes (see Fig. 
[TJ which captures the flash event we have modelled). This as- 
sumption will enable us to use some models at more than one 
sampling point, since in contrast to, e.g., T e g, the stellar mass M+ 
is not a very critical parameter in the wind and CSE models. For 
the luminosity and temperature this 'recycling' of models leads 
to deviations from the evolutionary track no larger than 1%. 
Exact matching of stellar parameters will not increase the accu- 
racy of the output, since small changes (< 1%) in the input pa- 
rameters have virtually no effect on the resulting wind and mass 
loss rate. We also assume a constant C/O ratio (where log(so) = 
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Fig. 1. Time evolution of luminosity L* (upper panel, full line), effective temperature r e (j (upper panel, dashed line), mass loss rate 
dM/dt assumed in the evolution model (lower panel, full line) and stellar mass M* (lower panel dashed line). The sampling points 
for the stellar parameters used as input to the dynamical atmosphere model are tick-marked in the lower part of each panel and 
the dotted line in the lower panel indicates the stellar mass that we adopt for the atmosphere models. The jump in T e ff seen near 
t — 1 177.58 Myr is a numerical artifact. 



-3.35) before and after the flash event, where C/O = 1 .85 before 
the flash, and C/O = 2.01 during and after the flash. In practice, 
this is really not an approximation, since on an evolutionary time 
scale the carbon abundance in the atmosphere makes a sudden 
jump from one basically constant level to another when a He- 
shell flash takes place. The numbers given above correspond to 
the evolutionary track we have computed. 

2.5. CSE Evolution 

Since the mass loss rate and wind velocity are changing with 
time, so are the boundary conditions for the modelling of the 
CSE around the star (the region where the detached shell is 
formed). The inner boundary conditions for such a model, i.e. 
wind velocity Vrj and density po at Rq (the radius where the mass 
loss is obtained) as functions of time, are defined by interpolat- 
ing the time sequence of atmosphere models computed as de- 
scribed above. Here, we are interested in whether these time 
variations will lead to the formation of a detached shell. The 
next step is therefore to compute the dynamical CSE evolution 
subject to such variable boundary conditions. 

Some combinations of stellar parameters and piston ampli- 
tudes Av p will result in models without any mass loss (see Table 
0. For such cases - if they occur in the post-flash phase - we 
assume a minimal mass loss rate of 10 M Q yr" 1 , which is moti- 
vated by observations (Olofs son et al. 1998llSchoier et al. 2001 , 
Scho ier et al. 20 05 ), and estimate the mimimum wind velocity 
to be ~ 1 km s . The circumstellar density at Rq can thus be es- 
timated using Eqn. (01. Another possibility would be to simply 
take the average density at Rq from the dynamical model with 



the the outer boundary radius Rq set closer (~ 1 .5/?+) to the pho- 
tosphere and assume a minimal mass loss rate, say 1O~ 8 M yr~' . 
Then the wind velocity can be estimated through Eq. (0. We pre- 
fer the first method, where we set vo = 1km s _1 , although both 
methods give more or less the same numbers. In the transition 
region between intense dust-driven mass loss and no mass loss, 
we set the rate to 1O~ 7 M yr~' and assume the wind velocity to 
be 3 km s and compute the density as described above. The 
exact numbers are not important in this case, since these mod- 
els represent a phase with essentially no stellar wind that occurs 
after the He-shell flash. 

We employ a simplified model of the CSE evolution where 
the equations of continuity and motion are given by 

% + V ■ (pv) = (4) 
at 

^-+V-(pvv) = -V/> gas , (5) 
of 

where p is the mass density, v is the velocity and P gas is the ther- 
mal gas pressure. These equations are solved using a simplified 
setup of the dynamic atmosphere code (radiative transfer and 
dust equations switched off), i.e., the numerical method is the 
same. We assume that the gravitational pull from the star does 
not affect the expansion of the CSE on the typical distance from 
the star considered here (R > 10 15 cm), since at this radius it 
is orders of magnitude weaker than in the atmospheric layers. 
The self-gravity of the CSE can also be considered negligible 
as the typical mass density is too low to produce any sigificant 
self-gravity effects. We also assume that the radiation field is 
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dominated by the central star, i.e. a point-source field falling off 
as ~ r~ 2 . Hence, at large distances from the central star, neither 
gravity, nor radiation, will be a force that really matters. This is 
of course a somewhat simplistic approach, but it is certainly a 
fair approximation of the physical situation in the present case. 
The pressure gradient, however, is kept in the equation of mo- 
tion, since it may play an important role in shock waves occuring 
at the interfaces of different wind phases. 

Regarding the temperature structure, we consider two limit 
cases. 

I: No radiative heating or cooling, the temperature structure is 
set by the expansion of the gas, leading to a fast decrease of 
temperature with distance from the star. Shocks are assumed 
to be adiabatic. 

II: A fixed, given temperature structure only reflecting the radi- 
ation field of the star, leading to a more gentle dropping of 
temperature with distance. Isothermal shocks (relative to the 
surrounding medium). 

More precisely, in Case I we are solving an energy equation 
of the form 



^ + V ■ (pe v) = -P gas V • v, 



(6) 



(where e is the specific internal energy and all other variables 
are denoted as before) together with the equation of motion and 
the continuity equation (see above). This means that the energy 
equation used for the CSE differs from the full RHD version 
of the atmosphere and wind acceleration region by the lacking 
source and sink terms for radiative heating and cooling (cf., e.g., 
Hofner et al.2003). 

In Case II, representing the other extreme, we assume that 
the temperature structure is fully determined by radiative equi- 
librium with the stellar radiation field. In this case, we have 



Ky(Jy - S V )dV = 0, 



(7) 



where k v is the monochromatic opacity, J v is the mean intensity 
and S v is the source function. Hence, the zeroth order moment 
equation of radiation, 



V ■H+p(k s J v -k s S) = 0, 



(8) 



where the frequency integrated opacities kj and ks are defined 
as in Hofner et al. (2003), implies that H in such a case must be 
divergence free, e.g. like a point source field. We thus estimate 
the radiative field from 



J v « W{r)B v {T±), 
where 



(9) 



W(r) = i 



is a geometrical dilution factor (cf. ILamers & Cassinelli 1999). 
If we assume radiative equilibrium for dust, and T gas » T^usn we 
end up with the relation 



where Qn(a, T) is the Planck mean extinction efficiency as a 
function of the grain radius a and the temperature. A simple esti- 
mate of the frequency dependency can be made from assuming 
a power-law relation Q(a, v) ~ v". In such a case, 



Tdust = W(r) 1/(4+B) TV, 
which for r » becomes 



Tdu 



R ^ v2/(4+n) 



(?) 



(ii) 



(12) 



Considering amorphous carbon grains, a reasonable choice for 
the power index is n « 1, i.e. a linear dependence on frequency. 
For the strictly linear case we get T^st ~ ^* r~ 2 ^ 5 , which is con- 
sistent with the temperature structure found by SS2000 using 
their DEXCEL-code. Hence, we have replaced the energy equa- 
tion in our CSE model by inserting this approximate temperature 
structure for the gas temperature 



?CSE 



-2/5 



(13) 



where Tq a 500 K and Rq — 7 ■ 10 cm. This approxmation is 
what we refer to as Case II, which is nearly the opposite to Case 
I, since the temperature structure is only reflecting the stellar 
radiation field and neither affected by shocks, nor by the expan- 
sion. It is reasonable to assume that T- d d < T ga% < Jdust, i-e., the 
true gas temperature structre will be somewhere in between the 
temperature structures implied by Case I and Case II. 

3. Results 

Stellar parameters and model results for the set of atmosphere 
models used to compute the mass loss history are presented in 
Table|2] A graphical display of the time evolution of the luminos- 
ity and the mass loss rate are shown in Fig. [2] where the expected 
correlation between these two quantities is evident. In Fig.[3]we 
compare the wind velocity and the mean degree of dust conden- 
sation - two quantities that show a correlation as well. In these 
plots, full lines represent the case where Av p = 4.0 km s _1 , while 
dashed lines correspond to Av p = 6.0 km s _1 . Increasing this pis- 
ton amplitude leads to a much stronger pre-fiash wind. The mass 
loss rate as well as the wind velocity changes significantly (as is 
evident from Fig. [2] and Fig. [3j. 

A close up of the evolution during the flash event (right pan- 
els in Fig.|2]and Fig. [3j reveals an anti-correlation between mass 
loss rate and wind velocity at the peak of the flash. This may 
seem curious, but can be explained by the fact that the relative 
abundance of dust is lower in the flash-peak model, i.e the model 
corresponding to the highest luminosity (see Tabled. A smaller 
abundance of dust (lower paust) w iH result in a less effective mo- 
mentum transfer between radiation and matter. Hence, we could 
expect a lower wind velocity, but not necesserily a lower mass 
loss rate, since the latter is proportional to both the matter density 
of the atmosphere and the wind velocity. In the particular case 
of the flash-peak model the gas density is higher, which leads to 
a higher mass loss rate although the wind is slower. The corre- 
lation between the mean degree of dust condensation f c and the 
mean wind velocity is seen in Fig. [3] 



TL (Ma, Tdust) = W(r) T* Q n (a, T+), 



(10) 



1 If the monochromatic Q is frequency-independent (gray case), the 
Planck mean Q?\ is independent of T and solving Eqn. d 1 0b for rd ust 
gives Tdust ~ T» r 1 ' 2 for r » R*. Opacities are, however, frequency 
dependent and Qpi(a, T) is therefore not a constant. 
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Table 2. Input parameters (L+, T e g, C/O, P, Av p ) and the resulting mean mass loss rate, mean velocity at the outer boundary and 
mean degree of dust condensation at the outer boundary. is the stellar radius of the hydrostatic initial model (obtained from L* 
and r e ff) and log(g) is the corresponding surface gravity. The dust-to-gas mass ratio pdust/Pgas is calculated from / c as described in 
Hofner & Dorfi (1997). All models have a mass of M* = 1.2M and log(eo) = -3.35. The mean values of the wind properties in 
this table are all taken over at least 100 pulsation periods, which is sufficient to ensure a small error due intermittent variations of 
these quantities. 



Model U r cff R, log(g) C/O P Av p (M) <v cxt ) (f c > flta 

[L G ] [K] [R e ] [days] [kms- 1 ] [M G yr 1 ] [kmr 1 ] [10" 4 ] 
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3.9 
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19 
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18 
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12000 
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19 


0.32 


18 



The fact that the mass loss rate, wind velocity and f c are 
strongly dependent on the piston amplitude Av p means that the 
predicted wind and mass loss cannot be very precisely con- 
strained. There is no way of knowing if this amplitude can be 
regarded as constant over time without a complete physical de- 
scription of the pulsation mechanism. This is certainly an impor- 
tant remark in the present case, since it appears quite likely that 
the pulsation amplitude could change as the star goes through a 
He-shell flash. Furthermore, the difference in kinetic energy in- 
jection by the pulsations depending on the choice of Av p may 
have significant effects on the amount of material being trans- 
ferred out to the dust-formation region, and thus also for the 
degree of dust condensation and consequently for the momen- 
tum transfer and wind velocity. For a further discussion on con- 
straints on Av p , see Section l4~T1 



3.1. Time Evolution of the Atmosphere 

Since the wind properties may show short-term, intermittent 
variations over time, the model output is time-averaged. The 
models usually reach a statistical steady state after ~ 50 pul- 
sation periods. An interval of ~ 100 periods is sufficient to ob- 
tain reliable mean values. A typical wind model is evolved over 
150 - 200 periods, which is thus sufficient for the output to be 
statistically stable. 



There is a complex interplay between density, dust forma- 
tion, radiation pressure and gravity (escape velocity) that makes 
the variations in mass loss for different sets of stellar parame- 
ters highly nonlinear. This interplay is why simple parameteri- 
sations, where the mass loss rate is essentially just a function of 
the luminosity (e.g. Blocker 19951 IWachter et al. 20021 1. cannot 
provide a satisfactory prescription for mass loss throughout the 
evolution of a star on the AGB. The well-known (but unfortu- 
nately often abused) Reimers formula (Reimers 1975) provides 
a mass loss prediction that cannot be made in agreement with 
the resulting mass loss evolution from our models. Regardless 
of the choice of the parameter r/R (in Fig. |4] the case 77R = 1 
is shown) the mass loss evolution is not as strongly peaked at 
the time of the He-flash as the evolution implied by our models. 
However, Reimers formula was intended as a description of the 
mass loss from red giants, which does not include TP-AGB stars, 
and there are no justifications for extending the use of Reimers 
formula to this regime. Its linear dependence on the luminosity 
L + means that it cannot at the same time reproduce the mass loss 
associated with the flash peak. A comparison of our results with 
existing mass loss formulae is presented Fig. |4] in which it be- 
comes clear that none of the simple descriptions considered here 
agrees with our detailed model, nor are they consistent with each 
other. It is particularly interesting to notice that the description 
by Wachter et al. (2002) yields a significantly higher over-all 
mass loss rate compared with other prescriptions, as well as our 
wind model. The reason for this, as far as we can tell, is two-fold: 
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Wachter et al. use grey radiative transfer (with lower gas opaci- 
ties) which is known to give higher mass loss rates and they de- 
rive their formula in principle only from models with saturated 
dust-formation, i.e. models in the regime where the momentum 
transfer is most efficient. 

In both Case I and Case II, this thin shell structure survives 
over the whole simulated CSE evolution. In Case II the thin shell 
structure remains quite intact as the shell expands - the only im- 
portant change over time is a slight broadening of the density 
peak likely due to internal pressure. The evolution of Case I is 
somewhat different. Here a very pointed shell structure evolves 
shortly after flash event, but as the shell expands two shock 
fronts are formed, creating a density bump (or plateau) altough 
the initial density peak remains. The density peak is located at 
the contact discontinuity between the slow (pre-flash) and fast 
wind phases, where the gas temperature reaches a minimum. 



3.2. Formation of a Detached Shell 

The results of the hydrodynamic CSE evolution modelling con- 
firm the "modified He-shell flash scenario" presented by SS2000 
in the case of Av p = 4.0 km s _1 . Using the mass loss evolu- 
tion we have obtained from our detailed modelling, we receive 
a density structure that is qualitatively consistent with the re- 
sults of SS2000, but in general results in a much thinner shell. 
Case I also shows a time variation of the relative thickness not 
present in the results of SS2000 (see Fig. [9] Fig. QT|in this pa- 
per and Fig. 12 in SS2000). The reason for this time variation of 
the thickness is probably the different evolution of the velocity 
profile shown in our results, shortly after the He-shell flash. As 
the CSE evolves, the velocity profile becomes very similar to the 
result from the "modified He-shell flash scenario" presented in 
SS2000. 

We have also calculated the FWHM of the quantity, dm/dr - 
4n r 2 p, which we will refer to as the thickness of the shell. In 
Fig. (0 we show the time evolution of the thickness of the shell 
for Case I and II. As is evident from the plots, including a rea- 
sonable temperature structure has really a dramatic effect on the 
CSE evolution. Beyond R = 10 17 cm (roughly 2000 years af- 
ter the flash event) the relative thickness AR S /R S of the shell 
appears constant in Case II, while in Case I we see more of 
a AR S /R S ~ l/R s behaviour. The typical relative thickness in 
Case II is AR S /R S ~ 0.01, which is less than what has been de- 
rived from observations (e.g. Olo fsson et al. 20001 . It should be 
pointed out, however, that these observations of detached shells 
we refer to are CO detections. Gonzalez-Delgado et al. (2001) 
have detected shell structures around AGB stars in the optical, 
which have AR S /R S ~ 0.1. This is an order of magnitude thicker 
than our model result and may indicate a problem with the as- 
sumptions made in Case II. In Case I the thickness is measured 
on the very sharp and pointed density peak (making the shell ap- 
pear extremely thin), which means that the numbers we derive 
are only of interest for comparison with Case II. 

The differences ocurring from the use of boundary condi- 
tions derived from atmosphere models with different piston am- 
plitudes, are quite dramatic. With Av p = 4.0 km s the CSE 
model provides a pronounced shell structure, while the Av p = 
6.0 km s case does not lead to the formation of any shell struc- 
ture at all (see Fig.[10]and Fig.fTTI). Since the Av p = 2.0 km s _I 
case was directly excluded (no wind, except at the flash peak) 
we are left with only one case that can produce a detached shell. 
In this case the density contrast at the shell radius is roughly two 
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Fig. 5. Timeseries of the wind velocity and mass loss rate at 
Ro showing the short-term variations before the He-flash (model 
180t28ml2cl85). 



orders of magnitude relative to the surrounding medium, which 
should definitely be a detectable feature. 

In order to determine if there are any effects on the shell 
formation from the short term variations of the wind during 
the flash (see Fig. [3j, we have integrated the mass loss dur- 
ing the flash event and created a representative pulse function 
(same integrated mass loss) and replaced the time dependent in- 
ner boundary condition with this simplified function (see Fig. 
|8}. The wind velocity is turned into a pulse function in the same 
fashion. Comparing the resulting CSE evolution with "case B" 
and "case C" in SS2000, it is clear that our CSE model produces 
a thinner and less massive shell, but still with an equally large 
density contrast. It is also evident that the small features (e.g. the 
dip at the flash peak) in the wind evolution found in our detailed 
modelling, have very little effect on the properties and evolution 
of the CSE and the detached shell (compare Fig.fTOland Fig.fTTI 
with Fig.fHland Fig.fBl). 

For the sake of completeness, we have tested how the 
mass loss history suggested by observations of TT Cygni 
dOlofsson etal. 19981 see also Fig. 1 and Fig. 4 in SS2000) af- 
fects the formation of a detached shell, using a boundary for 
the wind velocity which is equivalent to "case B" and "case C" 
in SS2000. As we show in Appendix lAl the results we get for 
Case II are fully consistent with the results of SS2000 using their 
"NEBEL" code. 

These tests elucidate the importance of modelling the wind 
evolution correctly. If the detailed inner boundary conditions de- 
rived from wind models are replaced by a combination of sim- 
ple pulse functions (conserving differences in mass loss rate and 
wind velocity, before, during and after the flash), the CSE evolu- 
tion is not much affected. If, on the other hand, the mass loss and 
velocity evolution differs as in "case B" and "case C" in SS2000, 
the CSE evolution will be different and results in a much thicker 
shell structure (e.g., compare Fig.[TT]with Fig. lA.2l and Fig. IA.4b . 
Thus, it is very important to use consistent combinations of mass 
loss rates and wind velocities when modelling the evolution of 
the CSE. 
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Fig. 2. Upper left panel: the external luminosity for Av p 
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time. Lower left panel: the mass loss rate for Av p = 4.0 kms -1 (full line) and Av p = 6.0 kms -1 (dashed line) as a function of time. 
The right panels shows a close ups (the interval between the vertical dotted lines) demonstrating the effects of the thermal pulse. 
The vertical dashed line marks the luminosity peak of the He-flash. 
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Fig. 6. Timeseries of the wind velocity and mass loss rate at Ro 
showing short-term variations during the peak of the He-flash 
(model 112t265ml2c201). 
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Fig. 7. Timeseries of the wind velocity and mass loss rate 
at /?o showing short-term variations after the He-flash (model 
180t28ml2c201). 



4. Discussion 

4.1. Mass Loss Rates 

The pulsational instability in the stellar interior provides a mech- 
anism for levitation of material into the outer atmosphere, i.e. 
the injection of kinetic energy from the pulsation related shock 



fronts as they propagate outwards may be sufficient to push ma- 
terial out to a radius where dust can begin to form. It is this phe- 
nomenon that enables the large mass loss rates associated with 
TP-AGB stars. 

As described in Section [2~2l the pulsation is simulated by ap- 
plying a piston boundary condition. Increasing the piston ampli- 
tude in the models from Av p = 2.0kms~ 1 to Av p = 6.0kms _1 
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Fig. 3. Upper left panel: the wind velocity for Av p = 4.0 kms -1 (full line) and Av p = 6.0 kms _I (dashed line) as a function of time. 
Lower left panel: the mean degree of dust condensation for Av p = 4.0 kms 1 (full line) and Av p = 6.0 kms _I (dashed line) as a 
function of time. The right panels shows a close ups (the interval between the vertical dotted lines) demonstrating the effects of the 
thermal pulse. The vertical dashed line marks the luminosity peak of the He-flash. 



can raise the density in the dust forming region by as much as 
one order of magnitude. Consequently, the amount of dust that 
forms may increase significantly as well and in turn lead to a 
more efficient momentum transfer and thus greater wind veloc- 
ities. One can in princple fine tune Av p in order to get the wind 
properties that lead to shell formation for a specific set of stel- 
lar parameters. However, even if it is difficult to put strong con- 
straints on the mass loss evolution without a reliable model for 
the pulsation mechanism, we can still limit the range of piston 
amplitudes that are reasonable in the present context. As we have 
seen, too small amplitudes do not lead to wind formation (ex- 
cept at the flash peak) and too large amplitudes will not provide 
the pronounced variations of wind properties needed to obtain 
the interacting wind phases that are crucial for shell formation. 
Hence, the possibility to choose the piston amplitude is restricted 
and we have found that the amplitude implied by previous com- 
parsions with observations (Av p = 4.0 km s~') is a good choice 
(e.g. Nowotny et al. 2005b[ l. Furthermore, models with this am- 
plitude give mass loss rates and wind velocities that matches 
observations of TP-AGB stars with known detached shells very 
well (e.g. Schoier et al. 2005). From recent observations of all 
objects with known detached molecular shells, Schoier et al. 
(2005) derived a post-flash mass loss rate (which at the same 
time gives a hint about the pre-flash rate) of a few times 1O~ 7 M 
yr~\ and ~ 1O~ 5 M yr~' associated with a flash event, which is 
almost exactly what we get (see Fig. [2j. The picture one gets is 
rather consistent - the pulsations are nicely modelled by a piston 
boundary condition with Av p ~ 4.0 km s~' even if the physical 
mechanism is not modelled. We have thus not considered tempo- 
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Fig. 8. Plots of the pulse functions used to replace the time de- 
pendent inner boundary condition. 



ral variations of the pulsational amplitude during the flash event, 
even if they may occur. 



Lars Mattsson et al.: Mass Loss Evolution and the Formation of Detached Shells around TP-AGB Stars 



11 




Fig. 4. Comparison with some existing mass loss formulae. The plot shows the mass loss evolution during the considered He-shell 
flash adopting mass loss rates from our atmosphere model (for two different values of Av p ), along with a simple analytic model 
dLamers & Cassinelli 19991 ), the empirical Reimers (1975) formula and the mass loss prescriptions derived from grids of numerical 
models by Wachter et al. (2002) and Blocker (1995). Due to the large differences we have plotted the logarithm of the mass loss 
rate. The right panel shows a blow-up of the interval marked with vertical dashed lines in the left panel. Note as well, that for the 
evolution computed in this paper, the values at the mass loss minimum are ad hoc values. 



The three types of models that we have combined provide a 
rather consistent connection between a He-shell flash and de- 
tached shells with essentially one adjustable parameter, Av p . 
However, there is one possible inconsistency that we have not 
mentioned so far. The evolutionary tracks for the present stel- 
lar model have been calculated using the mass loss prescription 
provided by Blocker (1995). As shown in Fig. [4] the mass loss 
rates predicted by our model differ from most simple mass loss 
prescriptions (formulae). A significantly different mass loss evo- 
lution will alter the late stages of stellar evolution (as pointed out 
by IBlocker 19951 1 and thus provide a different set of stellar pa- 
rameters as input to our model of the wind evolution. But this is 
most likely a minor problemns in terms of consistency, since the 
time interval we model is short on an evolutionary time scale and 
each one in the series of thermal pulses corresponds to slightly 
different sets of parameters, regardless of how the mass loss af- 
fects the stellar evolution. Moreover, the flash event itself, is not 
regulated by mass loss but by nuclear reaction rates. In our mod- 
elling we have picked one typical flash event that looked suitable 
in terms of getting a reasonable pre-flash mass loss rate, judging 
from previous experience of mass loss modelling. 

4.2. CSE Evolution and the Formation of a Detached Shell 

Let us now discuss how a detached shell may, or may not, be 
formed in the two limit cases of our CSE model (see Section 
\2.5\ due to the wind evolution obtained from our wind models. 

In Case I, where the global temperature structure is set by 
the expansion of the gas, the shock is basically adiabatic. The 
increase of the pressure leads to high temperatures (typically a 
few 10 3 K) in the shocked regions, which cannot be radiated 
away (since radiation is not included) and thus enforces a change 



in the gas motion that manifests itself by the build up of a shock 
wave in front of the ejected shell. This is why the shell evolves 
into a density bump with a density peak on top of it. There are 
several arguments against this kind of "adiabatic" evolution of 
the shell, where the most important one is that any shock front of 
this type must lose at least some energy through radiation. The 
high temperature we obtain inside the shocked regions as well 
as the low temperatures in the stationary pre-flash wind, indicate 
that Case I is quite unrealistic. 

Case II, where we have assumed a fixed temperature struc- 
ture, corresponds to having an isothermal shock front relative 
to its surroundings. This is in principle the opposite extreme 
of Case I, since the expanding shell is forced to have the same 
temperature as the surrounding medium, meaning that all excess 
energy is transferred away from the shock. Such an effective- 
ness in radiative cooling is quite unlikely, but Case II is proba- 
bly still closer to reality than Case I - especially regarding the 
global temperature structure. To make a truly realistic model 
including all possible energy sinks and sources in this case, it 
is probably necessary to include a full RHD description (with 
dust) of the problem and maybe even go from ID to 3D in or- 
der to correctly treat turbulence and dynamical instabilities (see 
Myasnikov et al. 2000) . 

The reverse shock and the forward shock in front of the shell 
are effects due to the quite large velocity jumps that build up on 
each side of the shell as it propagates outward. This feature is 
most prominent in Case I, due to the overestimation of compres- 
sional heating and the resulting thermal pressure gradients. 

As we pointed out in Section [3~2l changing the mass loss and 
wind velocity evolution into simple step functions (where the in- 
tegrals are preserved) does not notably affect the resulting shell 
structure. The reason is that the formation of a thin shell happens 
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Fig. 9. Thickness of the detached shells formed in Case I & II. The left panel shows Case I and the right panel shows Case II for 
Av p = 4.0 kms -1 . The dashed lines indicate fits to AR S /R S cc 1/R S . Each black dot represents a specific time in the evolution of the 
detached shell in accordance with all following plots of the velocity and density evolution. The dotted line marks the lower limit for 
the typical shell thickness derived by Olofsson et al. (2000) from CO observations. 



very much through wind interaction, and since we do not change 
the wind velocity of the pre-flash phase or at the mass loss max- 
imum, the wind interaction takes place in a similar fashion. 

We return now to the importance of having a correct wind 
model as the basis for the boundary variations in the CSE model. 
In SS2000 the inner boundary of the atmosphere model is pre- 
scribed such that a mass loss rate is assumed and the computa- 
tions start beyond the condensation radius where the wind veloc- 
ity is set to some initial value (see Steffen et al. 1998 for details). 
The velocity field is thus not derived from a fully self-consistent 
description of the stellar wind formation (although the equation 
of motion is solved, subject to the given boundary conditions). 
Moreover, the mass loss rate is not computed, but rather chosen 
according to a semi-empirical mass loss formula (Blo cker 1995b . 
For the CSE evolution to be correcly modelled, however, we 
must be certain that the variations of both the mass loss rate and 
the wind velocity are essentially correct and consistent, since the 
interaction of wind phases is perhaps the most important ingre- 
dient. 

Our detailed wind models (for Av p = 4.0 km s~') show that 
the wind velocity evolves from a slow phase (~ 10 km s -1 ) be- 
fore the He-shell flash, into a rather fast wind (~ 20 km s _1 ) dur- 
ing He flash event and then back to a slower wind phase again 
(see Fig.|2|. This is largely the same type of evolution used in the 
"modified He-shell flash scenario" in SS2000, although we have 
obtained it using a more sophisticated model, which provides 
both mass lass rate and wind velocity as functions of stellar pa- 
rameters. 

In the Av p = 6.0 km s _1 case we obtain a very different evo- 
lution compared with the Av p = 4.0 km s _I case, also qualita- 
tively speaking. The pre-flash mass loss rate is almost an order 
of magnitude higher and the wind velocity is a factor of two 
higher, providing almost no variations in the wind velocity be- 
fore and during the flash. The latter is the main reason why the 
Av'p = 6.0 km s case does not render a shell structure. 

How reliable are these results then? One may ask whether, in 
some part, the pre-flash wind could be different if drift/coupling 
between the gas and dust phases was to be introduced. The ques- 
tion is very relevant since these effects are treated in SS2000. 
At the flash peak this is not a problem, since high mass loss 
rate results in a higher CSE density which makes the effects of 



drift much less significant. We are likely to overestimate the pre- 
flash wind velocity due to the phase coupling. Hence, the pre- 
flash wind velocity would probably be even lower if drift was to 
be included, making the velocity jump (see Fig. [3]l even larger. 
Therefore, the inclusion of drift would probably not change the 
result qualitatively. 



The thickness AR S of the shell is hard to evaluate in Case 
I, since the density peak has almost no spatial extension com- 
pared with the scales considered here, and the bump is probably 
a somewhat artificial structure created by a dramatic pressure in- 
crease at the shock fronts since no energy is radiated away. It is 
therefore difficult to interpret the meaning of the thickness as de- 
fined by the FWHM. In Case II the characteristic absolute thick- 
ness of the shell increases almost linearly with the shell radius 
(which in turn grows close to linearly with time). This results 
in a constant AR s /7? s -ratio, in agreement with the result found 
by SS2000 in the "modified He-shell flash scenario". However, 
where SS2000 find a ratio ARJR S w 0.04, we find this number 
to be roughly 0.01, as we have pointed out previously in Section 
13.21 If the faster wind associated with the "mass loss eruption" 
runs into a less dense medium, as in our model, it means that 
less matter is swept up and a thinner shell may thus be formed 
out of this interaction. But as far as we can tell, the thinner shells 
can not be explained only by the weaker pre-flash wind that we 
obtain, nor can it be explain by any significant difference in the 
duration of the mass loss eruption. A possible, additional expla- 
nation might be the fact that the shocks are isothermal with their 
surroundings in our model, and therefore the shell structure is 
not substantially broadened by thermal pressure. But the condi- 
tions in the SS2000 model are also close to isothermal, which 
makes this explanation seem less likely. A perhaps more likely 
explanation arises from the fact that the wind velocity jump from 
pre-flash to the flash-peak is an important parameter governing 
the formation of a detached shell. In the present work this jump 
is roughly a factor of 2 larger than in SS2000, which may be a 
viable explanation. In general, however, our results confirm the 
connection between stellar evolution (the He-shell flash) and the 
formation of detached geometrically thin shells around TP-AGB 
stars. 
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Fig. 10. Nine instants in the evolution of the CSE. Velocity and density structure of the CSE in Case I. Upper panels: the velocity 
field for Av p = 4.0 kms -1 (left) and Av p = 6.0 kms -1 (right). Lower panels: the logarithmic density for Av p = 4.0 kms _I (left) and 
Av p = 6.0 kms -1 (right). 
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Fig. 11. Same as Fig. [TO] but for Case II. 



5. Summary and Conclusions 

We have used the RHD code for dynamic atmospheres by Hofner 
et al. (2003) in combination with a stellar evolutionary track of 
a M = 2M Q , Z = 0.01 star in order to calculate the evolution 
of mass loss and wind velocity during a He-shell flash. The ob- 
served wind and mass loss properties associated with detached 
shells are rather well reproduced by our models of TP-AGB star 
atmospheres using stellar parameters sampled from this evolu- 
tionary track. None of the other existing mass loss prescriptions 
for AGB stars we have considered give, quantitatively speaking, 
mass loss rates similar to those which we compute. However, 
constraints on the mass loss evolution derived from observations 
(Olof sson et al. 2000 Schoie r~et al. 20051 appear to agree better 
with our results than with these mass loss formulae. Our numer- 
ical model may therefore be considered as a reasonably realis- 
tic description of the connection between stellar parameters and 
mass loss rates for TP-AGB stars. 

We find that both mass loss rate and wind velocity increase 
during the flash event. This is in qualitative agreement with the 
"modified He-shell flash scenario" proposed in SS2000 and we 
thus conclude that geometrically thin detached shells are neither 
likely to be formed by simple two-wind interaction nor out of 
a "mass loss eruption" (with no variations in the wind velocity) 
alone. In fact, the formation of a slow and a fast wind combined 
with an eruptive mass loss associated with the fast wind turns out 



to be very critical. To use consistent combinations of mass loss 
rates and wind velocities is therefore most important in order to 
actually investigate a connection between He-shell flashes and 
detached shells. The amplitude of the internal pulsations of the 
star has a great effect on the wind velocities before and after the 
flash and thus on the strength of the wind interaction, in the case 
studied here. 

Our simplified dynamical model of CSE evolution agrees 
qualitatively with what is expected from observations. We con- 
sidered two basically complementary approximations for radia- 
tion effects. First we considered the case with no radiation effects 
at all (Case I) and then a second case where we used a simple 
analytic model for the temperature structure based on radiative 
equilibrium for the dust (Case II). In both cases we have assumed 
that all forces, except thermal pressure, are negligible. Case I and 
Case II agree qualitatively with each other, in the sense that they 
both produce a geometrically thin detached shell. However, only 
in Case II the relative thickness of the shell seems marginally 
consistent with CO observations. The thin shell structure seems 
to be well maintained throughout the evolutionary time span of 
the model, but the thickness of the shell grows linearely as it 
expands, which is most likely an effect of the internal pressure. 
Case I initially forms a shell structure with a much more pointed 
appearence, which in fact is the contact discontinuity at the inter- 
face between the winds. The pointed structure seems stable, but 
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Fig. 12. Nine instants in the evolution of the CSE. Velocity and 
density structure of the CSE using the step function boundary 
condition (see Section |3~2l i in Case I (Av p = 4.0 kms _I ). 
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Fig. 13. Same as Fig.Ql] but for Case II. 



as the expansion continues a reverse shock propagating inwards 
is created, as well as an outgoing shock front which sweeps up 
a large portion of the matter ejected by the pre-flash wind. As a 
consequence of this, the shell structure evolves into a wide "den- 
sity bump" with an unreasonably thin shell structure on top of it. 
From this we conclude that the temperature structure is one of 
the key components for the shell formation. In addition, simple 
ID models may be inadequate since a full 3D simulation may 
result in much more complex structures that, in turn, may have 
significant effects on the formation and evolution of the shell. 

Further work is obviously needed, but so far we can conclude 
that the He-shell flash and the associated "mass loss eruption" in 
combination with wind interaction gives a satisfying and consis- 
tent description of the formation of detached shells around TP- 
AGB stars. Our detailed modelling of the mass loss evolution 
shows that a very short, intense period of mass loss can be asso- 
ciated with a He-shell flash. The "mass loss eruption" is followed 
by a period of no or very little mass loss, before a steady wind 
phase is once again established. This "eruptive" mass loss evolu- 
tion does not seem to lend much support for the simple two-wind 
interaction scenario where a constant faster wind, with a higher 
mass loss rate associated with it, runs into a slower wind. Since 
a combination of a "mass loss eruption" and interacting wind 
phases is consistent with both stellar evolution and the observed 
properties of detached shells, we find this scenario to be the most 
likely. 



This paper has presented a first attempt to construct a self- 
consistent model of the formation of detached shells. It is, how- 
ever, important to consider the CSE evolution in more detail than 
we have done here. A full RHD treatment in 3D would proba- 
bly show additional features (e.g. instabilities) that might help to 
constrain the mass loss history in more detail. Our intention is to 
return to this issue in a future publication. 
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Fig. A.l. Nine instants in the evolution of the CSE. Upper panel: 
velocity as function of radius for Case I with the same boundary 
condition as "case B" in SS2000. Lower panel: density as func- 
tion of radius for Case I with the same boundary condition as 
"case B" in SS2000. 
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Fig. A.2. Same as Fig. IA.ll but for Case II. 
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Fig. A.3. Nine instants in the evolution of the CSE. Upper panel: 
velocity as function of radius for Case I with the same boundary 
condition as "case C" in SS2000. Lower panel: density as func- 
tion of radius for Case I with the same boundary condition as 
"case C" in SS2000. 
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Fig. A.4. Same as Fig. [A3] but for Case II. 



Appendix A: Comparison with case "B" and "C" in 
SS2000 

This appendix presents test cases using basically the same 
boundary variations as in SS2000. Their "NEBEL" code con- 
tains a detailed treatment of radiation and cooling, compared 
with our fixed temperature structure in Case II. The "DEXCEL" 
runs have no directly corresponding case among the ones we 
have tried here. But since the "DEXCEL" code is known to suf- 
fer from numerical diffusion (which was pointed out in SS2000) 
any direct comparison would be difficult. Nonetheless, we think 
it is important to show that our code (Case II) produces results 
that are comparable to the "NEBEL" runs presented in SS2000 
if the input is similar. 

Using a simple mass loss evolution of the kind used in 
SS2000 (i.e. step functions, as in Section 13.21 ) in our Case II, 
which is quite similar to the "NEBEL" models of SS2000, we 
are able to basically reproduce their results, i.e. the wind veloc- 
ity field and the mass density of the CSE (Fig. lA.2l and Fig |A.4b 
evolve in ways that are very similar to their "case B" and "case 
C". In Case I we obtain a result that shows some resemblance to 
the result from the "DEXCEL" model of SS2000, but we see no 
obvious signs of numerical diffusion. The corresponding plots 
for Case I (Fig. IA.1 l and Fig [A3) are shown for comparison. 



